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We report on our investigation of the crystal structure of arsenic under compression, focusing 
primarily on the pressure-induced A7 — > simple cubic (sc) phase transition. The two-atom rhombo- 
hedral unit cell is subjected to pressures ranging from GPa to 200 GPa; for each given pressure, 
cell lengths and angles, as well as atomic positions, are allowed to vary until the fully relaxed struc- 
ture is obtained. We find that the nearest and next-nearest neighbor distances give the clearest 
indication of the occurrence of a structural phase transition. Calculations are performed using the 
local density approximation (LDA) and the PBE and PW91 generalized gradient approximations 
(GGA-PBE and GGA-PW91) for the exchange-correlation functional. The A7 sc transition is 
found to occur at 21±1 GPa in the LDA, at 28±1 GPa in the GGA-PBE and at 29±1 GPa in the 
GGA-PW91; no volume discontinuity is observed across the transition in any of the three cases. We 
use fc-point grids as dense as 66x66x66 to enable us to present reliably converged results for the 
A7 — > sc transition of arsenic. 



I. INTRODUCTION 

It is currently de rigueur to examine, both experimen- 
tally and theoretically, the structure of materials sub- 
jected to extremely high pressures li The group- V semi- 
metals, which include arsenic, are of particular interest 
as they all exhibit a distinctive low-symmetry structure 
when uncompressed. With pressure, these elements ex- 
perience transitions into structures of higher symmetry, 
with unusual intermediate phases — incommensurate or 
host-guest structures — occurring along the way. 

The computational study of the high-pressure behav- 
ior of materials necessitates the ability to properly ex- 
amine displacive phase transitions from one structure to 
another. Some structural phase transitions are straight- 
forward enough that the transition pressure can be deter- 
mined by looking at where the enthalpy-pressure curves 
of the two structures cross, or by using the technique 
of the "common tangent" on the energy-volume curves 
of the two structures. However, the transition pressure 
cannot be determined in this way if the enthalpy-pressure 
(or energy- volume) curves of the two structures merge, as 
in the case of the A7 — > sc transition of arsenic (and for 
the same transition in the other group-V semi-metals). 
So if the phase transition between the two structures is a 
smooth one — such that the energy differences between 
the two structures in the region of the transition are ex- 
tremely small — then how should such a transition be 
studied? High levels of convergence of the calculation of 
the quantities of interest are crucial. 

In the past, the A7 —^ sc transition pressure of ar- 
senic has been the subject of experimental as well as 
theoretical dispute. Experimentally, there has been a 
long-standing question as to whether results obtained 
by Beister— for this transition pressure are correct over 
those of Kikegawa and Iwasaki^ — our results support 



the findings of the former. Existing theoretical stud- 
ies of arsenic yield a wide range of possible values for 
the transition pressure. We believe that a great part of 
the reason for this spread of values is inadequate fc-point 
sampling of the Brillouin zone. The Fermi surface of ar- 
senic is extremely complex: using values reported by Lin 
and Falico\ii for the cross-sectional area of a neck (the 
finest feature of the hole Fermi surface), we estimate that 
to sample the Fermi surface of arsenic at ambient pres- 
sures such that all of the features can be resolved — and 
are subsequently correctly weighted — a grid of at least 
140x140x140 fc-points would be required. In practice, 
it is not necessary to use grids as dense as this but one 
would not know this without proper and thorough con- 
vergence testing with respect to fc-point grid size and 
amount of smearing used. As no such convergence test- 
ing of the A7 sc transition of arsenic appears in the 
literature, results of computational studies to date are 
unreliable and any agreement with experiment has been 
fortuitous. We undertake indepth studies of the conver- 
gence properties of this transition;^ which enable us to 
proffer reliably converged results. 

This investigation is not only timely, given the current 
levels of interest in the field of high-pressure research, 
but it is also especially relevant to the study of pressure- 
induced insulator to metal or semi-metal to metal phase 
transitions. In particular, we find that inspection of such 
transitions demands careful convergence of the quantities 
of interest with respect to fc-point sampling and smear- 
ing, yet such testing is not always performed. 

Using density functional theory (DFT) to study the 
high-pressure properties of materials furthermore necessi- 
tates an awareness of how different approximations to the 
exchange-correlation functional perform at such a task; 
to understand the differences encountered in the use of 
these functionals is in itself a motivation for undertaking 



theoretical investigations of materials at high pressures)^ 
We thus compare the performance of the LDA with that 
of two GGA functional, having achieved accurate results 
for each. 

The paper is organized as follows. Section II contains 
the background material that pertains to our study, in- 
cluding a description of the A7 phase, a review of the 
structural transformation that yields the sc phase when 
a sufRciently high external pressure is applied, and a sur- 
vey of the relevant findings published prior to this work. 
Section III deals with the computational details of our 
calculations. In Sec. IV, we present and discuss our re- 
sults, contrasting them against those that appear in the 
literature. We begin by examining, using the nearest and 
next-nearest neighbor distances, the overall picture of the 
pressure-dependence of the phases of arsenic for the LDA 
and the GGA cases; the pressure-dependence of the lat- 
tice parameters is also explored. Next we focus on the 
A7 — s- sc transition in particular, discussing it in more 
detail using energy-volume and pressure-volume curves. 
In Sec. V, we provide a brief discussion and conclusion. 



II. THE STRUCTURES OF ARSENIC 

At ambient pressures, arsenic is a covalently bonded 
compound existing in the rhombohcdral (A7) phase. 
The A7 phase, belonging to space group 166, is a 
low-symmetry, three-fold coordinated, layered structure 
(Fig. [T]). Nearest neighbors exist within a common layer; 
bonding between nearest neighbors gives the structure 
its buckled appearance. Adjacent layers, which contain 
next-nearest neighbors, are weakly bonded. The stack- 
ing of the layers occurs in the [111] direction, as does the 
buckling. 

The rhombohcdral primitive cell of arsenic contains 
two atoms and is described by the length of the prim- 
itive lattice vectors, a, the angle, a, between each set 
of primitive lattice vectors, and by the atomic positional 
parameter, z, which determines the positions of the two 
atoms along the cell's body diagonal. The two atoms are 
located at (Wyckoff) positions (z, z, z) and (— z, — z, —z)r 
along the primitive lattice vectors, where the atomic po- 
sitional parameter z is a fractional coordinate. Together, 
a and z characterize the degree of buckling experienced 
by the A7 structure. 

The A7 structure can be interpreted as a distortion 
of the six-fold coordinated sc structure^' '^"'^^ a diagram 
charting the evolution of the sc structure into the A7 
structure can be found in Ref. [3- The sc lattice can be 
viewed as being composed of two interpenetrating face- 
centered cubic (fee) lattices. The vertices of a sc cell 
along the [111] direction contain one atom from each of 
the fee lattices. A rhombohcdral cell is superimposed 
onto this in such a way that its origin lies central to, 
and contains, two of these alternate fee lattice atoms. Its 
body diagonal is along the [111] direction (the trigonal 
axis). Along this direction, the atomic positional param- 
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FIG. 1: (color online) The [111] direction is out of the page 
in the top portion of the figure; in the bottom portion of 
the figure it lies in the plane of the page and points upward, 
(a) Sizes and colors of atoms depicted carry no significance, 
but are employed merely to illustrate the ABC stacking of 
the A7 structure of arsenic, (b) The A7 structure consists of 
layers of 6-member buckled rings as displayed by the hexagon 
composed of green arsenic atoms. The buckled appearance 
of arsenic is due to the bonding of nearest neighbors that 
occurs within each of the layers; the layers themselves are 
only weakly bonded. Nearest and next-nearest neighbor dis- 
tances are illustrated using the red atoms and are measured 
with respect to the black atom. Nearest neighbor distances 
exist between red and black atoms within the same layer; 
next-nearest neighbor distances exist between red and black 
atoms from adjacent layers. The two-atom rhombohedral unit 
cell is also displayed. (This figure was created in part using 
AtomEye.-) 



eter z can be thought of as the ratio of the distance from 
the origin to the first atom, divided by the length of the 
body diagonal. A shear is then applied along the trigonal 
axis so that the angle a between sets of primitive lat- 
tice vectors decreases from its simple cubic value of 60°. 
Finally, the two atoms are displaced from their original 
positions and away from each other such that z decreases 
from its sc value of 0.25. The periodicity of the atoms in 
the [111] direction is now doubled; pairs of atoms in this 
direction repeat a pattern of being further apart and then 
nearer. In this way, the coordination number is reduced 
from six to three, as is consistent with the tendency of 
group-V semi-metals at normal pressures to form three 
bonds. 

Subjected to sufficiently high pressures, arsenic under- 
goes a transition from semi-metallic A7 to metallic sc. At 
lower pressures, the stability of the A7 structure against 
that of the sc structure can be understood in terms of a 
Peierls-type distortion associated with the displacement 
of the atoms along the [111] direction from their original 
positions within the simple cubic lattice. This displace- 
ment corresponds to the longitudinal acoustic phonon 
mode at the point R = 7r(l,l,l)/a on the boundary 
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of the simple cubic Brillouin zonei^^ The breaking of 
symmetry causes a lowering of the energy of the elec- 
tron states near the Fermi level, with the result that a 
gap is opened up, the band contribution to the total en- 
ergy is reduced, and the structure is stabilized against 
the distortion. As the pressure is increased, ionic effects 
become more significant, and the Ewald contribution to 
the total energy increases. Changes in the Ewald energy 
eventually dominate over changes in the band energy, 
the A7 structure is destabilized such that the Ewald en- 
ergy is lowered, and the sc structure is recovered f^SJ^ 
Within the framework of chemical bonding, at low pres- 
sures the nearly orthogonal bonds of the A7 structure 
indicate weak s-p hybridization, where valence electrons 
arc localized within saturated covalent bonds. The split- 
ting between s and p bands increases with pressure, until 
the s band is so deeply depressed with respect to the p 
band that unsaturated orthogonal bonds form, resulting 
in the metallic sc structureJi^^ 

It is sometimes convenient to use the nonprimitive 
hexagonal (trigonal) representation of the unit cell of ar- 
senic. The trigonal cell is described by the lattice vectors 
a hex and Chox, the trigonal axis which lies in the [111] 
direction and coincides with the body diagonal of the 
rhombohedral primitive cell, and by the same atomic po- 
sitional parameter z. The six atoms contained within the 
cell are located at Wyckoff positions ±(0, 0, 0)r/i<^ The 
lattice vector lengths a hex and Chcx are related to the 
length of the primitive lattice vector a and angle a as fol- 
lows: ahex = 2asin(a/2), and Chox = a{3[l-l-2cos(a)]}^/^. 
Within this representation, z can be thought of as locat- 
ing, as a fraction of Chcx, the first atom from the origin 
along this axis. 

From ambient conditions and with increasing pressure, 
the experimentally observed transition sequence for the 
phases of arsenic is A7 — > sc ^ As-III bcc, where bcc 
is the body-centered cubic phase and As-III is an in- 
commensurate structure. In the literature, the values 
reported for the pressure at which occurs the first struc- 
tural transition of arsenic, from semi-metallic rhombo- 
hedral to metallic simple cubic, arc spread out over a 
range of approximately 20 GPa. It has been theoreti- 
cally predicted that arsenic should undergo a structural 
phase transition at 18 GPa^^ early numerical studies 
of arsenic yielded for the transition pressure values of 
35 GPapi^ and 19 GPaf^i This last result corresponded to 
VjVo — 0.8 (and was indeed found again very recently 
in a study of the lattice dynamics of arsenic to be the 
volume ratio at which a transition occurs^), where V is 
the volume of the compressed cell and Vo is the uncom- 
pressed equilibrium volume. Investigating the pressure- 
induced superconductivity of arsenic, the experimental 
findings of Kawamura and Wittigi^ suggest a transition 
at 24 GPa, while those of Chen, et alt^ support a tran- 
sition pressure of 32 GPa, though they obtain 36 GPa 
via the theoretical component of their study. Perform- 
ing energy-dispersive and angle-dispcrsivc powder X-ray 
diffraction studies of arsenic up to 45 GPa, Kikegawa 



and Iwasaki^ determined the A7 sc transition to be 
of first order, occurring somewhere in the interval of 
31.4-36.6 GPa and resulting in a cell volume disconti- 
nuity of AV/Vt — 5%, where AV is the difference in 
the cell volume across transition, and Vt is the volume 
of the cell just prior to the transition. Beistcr, et al.^ 
conducting an angle-dispcrsivc powder X-ray diffraction 
study up to 33 GPa in addition to a Raman investigation, 
concluded that the transition occurs at 25 GPa, and that 
to within their experimental uncertainty no discernable 
change in volume is seen to occur: AV/Vt < 0.5%. More 
recently, theoretical investigations of the A7 sc transi- 
tion of arsenic yielded for the transition pressure a range 
of 25 ± 8 GPai^ in one study, and 28 GPai^ in another. 
A molecular dynamics simulation of the A7 sc tran- 
sition performed by Durandurdu^S yielded a transition 
pressure of 35 GPa and a volume discontinuity of 3.2%, 
and most recently, transition pressures of 22 GPa^i and 
26.3 GPa3^ have been found, with volume discontinuities 
of 0.8% and 0.4% respectively. 

Experiments aimed at investigating the higher- 
pressure phase transitions of arsenic were first carried 
out by Greene, et al.^ Energy-dispersive powder X-ray 
diffraction experiments on arsenic up to 122 GPa yielded 
sc — > As-III and As-III bcc transition pressures of 
48±11 GPa and 97±14 GPa, respectively. To within 
their experimental error of 1%, no measurable volume 
discontinuities for either phase transition were observed. 
These results were reproduced rather closely in a theo- 
retical study by Haussermann, et al.^ in which the two 
pressures were found to be 43 GPa and 97 GPa. Rein- 
terpration of diffraction patterns obtained in the experi- 
ments of Greene, et al. led Iwasaki^i (see also Ref. [251 ) to 
suggest a tetragonal Bi-III-type or Sb-II-type structure, 
regarded as distorted bcc, for the As-III phase. More 
recently McMahon, et al.^ using angle-dispersive X-ray 
diffraction techniques, determined both the Bi-III and 
Sb-II phases to be Ba-IV-type structures, consisting of a 
body-centered tetragonal (bet) host structure containing 
a bet guest component that is incommensurate with the 
host along the tetragonal c axis. They suggested that a 
modification of this structure in which the host compo- 
nent is bet and the guest is monoclinic would be appro- 
priate for As-III. Degtyareva, et al^ later found that 
the As-III structure closely resembles that of Sb-IV, and 
proposed that it is possible that arsenic may undergo an 
incommensurate-to-incommensurate transition between 
the As-III and bcc phases similar to the Sb-IV Sb-II 
transition observed in their study. 



III. COMPUTATIONAL DETAILS 

All calculations were performed using the CASTEP 
code.— For the exchange-correlation functional, we use 
the local density approximation (LDA)^ as parametrized 
by Perdew and Zunger^ the Perdew-Burke-Ernzerhof 
generalized gradient approximation (abbreviated as 
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GGA-PBE)2i and the Pcrdew-Wang generalized gradi- 
ent approximation (GGA-PW91)i^ 

The interaction between the core and valence states 
of arsenic is described using scalar-relativistic ultra-soft 
pseudopotentialsi^ For the calculations to proceed ef- 
ficiently, such pseudopotentials require the use of two 
FFT-grids: a standard grid, which is used to represent 
the pseudo wave- functions, and a fine grid, which is used 
to store the charge density, and thus the ultra-soft aug- 
mentation charge. The size of the standard grid is de- 
termined by the cut-off energy — the maximum energy 
accounted for in the plane-wave expansion of the wave- 
functions; the fine grid must be chosen dense enough such 
that the charge-density within the atomic cores is repro- 
duced with sufficient detail — this is a requirement to 
ensure well-converged forces and stresses. Convergence 
tests on our system at atmospheric pressure led us to 
choose a cut-off energy of 450 eV; we chose our fine grid 
to be twice as dense as our standard grid. With these 
choices, our energies, forces and stresses are converged 
to within 0.003 cV/atom, 0.001 cV/A and 0.01 GPa, re- 
spectively. 

We apply pressures between GPa and 200 GPa to 
the two-atom rhombohedral unit cell of arsenic, allow- 
ing cell lengths and angles in addition to atomic posi- 
tions to vary. For each applied pressure, the enthalpy 
is minimized such that the system is relaxed from the 
structure that exists at atmospheric pressure as deter- 
mined experimentally by Schiferli^ The geometry opti- 
mizations are made to continue until the differences in 
the forces and stresses between iterations are less than 
0.005 eV/A and 0.05 GPa, respectively. Proper con- 
vergence studies in electronic structure calculations are 
important and the lack of such testing can lead to severe 
inaccuracies 1^ Accordingly, we have thoroughly investi- 
gated the structural convergence of our system in the 
region of the A7 — *■ sc phase transition; we have exam- 
ined the convergence with respect to fc-point grid size 
and amount of cold-smearing^^ used — the results of 
this convergence testing can be found in the supplemen- 
tary material provided;^ In the vicinity of the A7 sc 
transition, the fc-point grids that have been studied are 
nxnxn, where n is 24, 25, 26, 33, 50 and 66 (all fc-point 
grids are unshifted, so that the gamma point has been in- 
cluded in the calculations), using cold-smearings of 0.1, 
0.2 and 0.5 eV in the LDA and 0.1 eV in the GGA-PBE. 
At pressures outside the region of this phase transition 
(and for all pressures studied in the case of the GGA- 
PW91); a fc-point grid of 33x33x33 is used, along with 
a cold-smearing of 0.1 eV. The symmetry of the initial A7 
structure is maintained throughout the relaxation, as a 
highly symmetric path is assumed for the A7 — > sc tran- 
sition, with atoms moving only along the [111] direction. 
Although it is irrelevant whether this assumption is valid 
for the As-III phase, because no attempt has been made 
here to properly model this incommensurate structure, 
it seems reasonable that imposing these constraints on 
the symmetry of the system will not affect the transi- 



tion pressures observed since the incommensurate phase 
is bounded by two phases of higher symmetry than the 
original A7 structure. 

Finally, all calculations are performed at T = OK, and 
the zero-point vibrational energy, as well as the spin-orbit 
interaction, are neglected. 



IV. RESULTS AND DISCUSSION 
A. Overview of the transitions of arsenic 

We first take an overview of the behavior of the near- 
est and next-nearest neighbor distances as the pressure 
is increased from GPa to 200 GPa. It has already 
been stated that for the A7 sc phase transition in 
particular, the behavior of these two quantities gives 
the clearest indication of when a structural phase tran- 
sition has occurred. Fig. [2] reveals the results we ob- 
tained for each of the exchange-correlation functionals 
studied, and includes separate close-ups of the A7 sc 
phase transition for the LDA and for the GGA. In- 
deed, it is quite clear when this transition occurs in all 
cases, as it happens when the nearest and next-nearest 
neighbor distances become equal. From the insets of 
Fig. m we can conclude that the A7 sc phase tran- 
sition of arsenic occurs at 21±1 GPa in the LDA, at 
28±1 GPa in the GGA-PBE and at 29±1 GPa in the 
GGA-PW91 (note that the difference in the transition 
pressure that arises from the use of the GGA-PW91 over 
the GGA-PBE is negligible as it is within the uncer- 
tainty of 1 GPa). This result is more consistent with the 
experimental findings of Beister, et al.^ who found the 
phase transition to occur at 25±1 GPa, than it is with 
those of Kikegawa and IwasakifS, who found it to occur 
in the range of 31.4-36.6 GPa; it furthermore disagrees 
with the outcome of simulations performed recently by 
Durandurdui^ 

From this same figure we note incidentally that the 
GGA nearest and next-nearest neighbor distances are 
always larger than those resulting from the LDA. This 
is due to the well-known fact that the LDA tends to 
overbind systems, whereas the GGA tends to underbind 
them<^ Thus cell volumes, cell parameters, and distances 
contained within the cell will tend to be lower for the 
LDA than for the GGA. 

Our unit cell contains two atoms and is periodically 
repeated. We cannot hope to say anything about the in- 
commensurate phase that occurs in the pressure-induced 
progression of phases of arsenic, and it is not our in- 
tention to do so in this paper; to acknowledge this 
point, we will use the abbreviation incomm. to rep- 
resent what we find between the sc and bcc phases, 
and subsequently write the progression of phases as 
A7 — > sc — s- incomm. — » bcc. We mention as an aside 
that we have chosen to work with the two-atom unit 
cell of arsenic in the interest of undertaking efficient cal- 
culations, but as we are guided by experiment and not 
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FIG. 2: (color online) Nearest neighbor (solid symbols) and next-nearest neighbor (open symbols) distances as a function of 
pressure for the LDA (black squares), for the GGA-PBE (red circles) and for the GGA-PW91 (blue triangles). The A7 — > sc, 
sc — > incomm. and incomm. bcc transitions are evident, where the insets emphasize the A7 — > sc transition in both the 
LDA and the GGA. (Calculations performed using a 33x33x33 fc-point grid and a cold-smearing of 0.1 eV.) 



primarily focused in this study on the incommensurate 
structure, a two-atom unit cell is sufficient for our needs. 
It is possible that were we to use unit cells containing ar- 
bitrary numbers of atoms we might discover other more 
complex or more closely-packed structures of arsenic at 
pressures higher than have yet been attained experimen- 
tally — this is an interesting point and although out- 
side the scope of this work, would be an excellent can- 
didate for the application of fully randomized structure 
predictioni^ 

Having stated that we do not expect to reproduce the 
incommensurate phase of arsenic, we can still use our 
system to approximate the pressure at which the system 
ceases to be sc, as well as to approximate the pressure at 
which it becomes bcc, and compare with the most recent 
experimental findings. Thus, referring again to Fig. [21 we 
see that in the LDA the sc —f incomm. phase transition 
appears to occur at 60-61 GPa and the incomm. — ^ bcc 
at 112-113 GPa. In the GGA-PBE it is 53-54 GPa and 
95-96 GPa, and in the GGA-PW91 it is 56-57 GPa and 
109-110 GPa. Although our findings for the sc in- 
comm. transition are about 10 GPa higher than those of 
Ref. [l^, and fall into the high-end of the range of values 
given in Ref. |23| . our results for the incomm. — > bcc tran- 
sition agree closely with those from both works. Along 
with the findings that are present in the literature, our 



results for the transition pressures of arsenic are summa- 
rized in Table m 

We break down the nearest and next-nearest neigh- 
bor distances into their constituent components in Fig. [31 
which illustrates the behavior of the lattice parameter a 
(top tier), the cell angle a (middle tier), and the atomic 
positional parameter z (bottom tier), for the LDA and 
the GGA-PBE. Each tier of Fig.[31is composed of a close- 
up of the A7 — > sc phase transition (on the left) along 
with the behavior of the quantity over the entire range of 
pressures studied (on the right). Respectively, the high 
symmetry values of a and z, as can be confirmed in this 
figure, are 60° and 0.25 in the simple cubic case, and 90° 
and 0.25 in the bcc case. The sc phase is of higher symme- 
try than the A7 phase, and because we have constrained 
the symmetry of the system, once the sc phase has been 
reached the atomic positional parameter z remains equal 
to 0.25. Thus, after the A7 sc phase transition has 
occurred, it is the angle a that determines the phase of 
the structure as the applied pressure is increased. 

The first point to observe upon inspection of Fig. [31 is 
that the A7 — > sc transition appears to be continuous 
with respect to all three quantities: a (top), a (mid- 
dle) and z (bottom). We see from the set of right hand 
panels in the figure that the sc — > incomm. and the in- 
comm. — » bcc transitions appear to be discontinuous with 
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FIG. 3: Lattice parameter a (top tier), cell angle a (middle tier) and atomic positional parameter z (bottom tier) as a function 
of pressure for the LDA (open circles) and for the GGA-PBE (filled squares). The left-hand panels emphasize the A7 — > sc 
transition, while in addition the sc — > incomm. and mcomm. bcc transitions are apparent in the top two tiers of the 
right-hand panels. (Calculations performed using a 33x33x33 fc-point grid and a cold-smearing of 0.1 eV.) 



respect to lattice parameter a and cell angle a, for both 
the LDA and the GGA; in particular, these discontinu- 
ities subsequently manifest themselves in the nearest and 
next-nearest neighbor distances at these two transitions 
as can be seen in Fig. and indicate when the cell an- 
gle has jumped away from 60° in the case of the former 
transition, and to 90° in the case of the latter. Before 
discussing the behavior with pressure of the cell param- 
eters further, we first inspect more closely the A7 sc 
transition. 



B. The A7 



sc transition of arsenic 



Examining the left hand panels of Fig. [3l we see 
that a and z reach their high-symmetry values at lower 
pressures in the case of the LDA than in the GGA; 
this is to be expected as we noted earlier that the 
LDA A7 sc transition pressure is lower than the 
GGA transition pressure. However, we also note that 
in both cases, a reaches 60° at a slightly lower pres- 
sure than z reaches 0.25. This result was also found by 
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TABLE I: Transition pressure Pt, reduced transition volume 
Vt/Vo, where Vt is the volume just prior to the transition and 
Vo is the uncompressed equilibrium volume, and percentage 
volume change across the A7 —> sc transition. Below are listed 
transition pressures for the higher pressure transitions of ar- 
senic. Parametrizations of the generalized gradient approx- 
imation to the exchange-correlation functional that appear 
below are PBB^ and PW91."^ (Abbreviations: PP - pseu- 
dopotential approach, LAPW - linear augmented plane wave 
approach.) 



A7 sc 


Pt (GPa) 


Vt/Vo 


AV/Vt (%) 


Theory: 








This work, LDA 


21 ± 1 






This work, PBE 


28 ± 1 
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Ref. 12 (PP, LDA) 


35 


0.72 




Ref. 11 (LAPW, LDA) 


19 


0.8 




Ref. 17 (PP, LDA) 


36 






Ref. 18 (PP, LDA) 


25 ±8 


0.78 


1 


Ref. W_ (PP, PW91) 


28 




0.8 


Ref. 20 (PP, LDA) 


35 




3.2 


Ref. 21 (PP, PBE) 


22 




0.8 


Ref. 15 (PP, PW91) 




0.8 
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26.3 




0.4 


Experiment: 








Ref. 


18 






Ref. 


24 






Ref._3 


31.4 - 36.6 


0.744 


5 


Ref. 2 


25 ± 1 


0.772 


0.5 


Ref.Tz 


32 






Pt (GPa) 


sc incomm. incomm. bcc 


Theory: 








This work, LDA 


60.5 ±0.5 




112.5 ±0.5 


This work, PBE 


53.5 ±0.5 




95.5 ±0.5 


This work, PW91 


56.5 ±0.5 




109.5 ±0.5 


Ref. IS (PP, PW91) 


43 




97 


Experiment: 








Ref. 23 


48± 11 




97± 14 









Seifert, et al^ in their study of the A7 sc transition 
of antimony. In fact, there is no reason to assume that 
a and z must reach their high-symmetry values at the 
same pressure — yet we cannot rule out the possibility 
that this pressure discrepancy may be an artifact result- 
ing from the use of smearing in our calculations. The 
difference in pressures between which a and z reach their 
high-symmetry values is slightly more pronounced for 
the LDA (approximately 3 GPa) than it is for the GGA 
(approximately 2 GPa for both the GGA-PBE and the 
GGA-PW91). The curves of the lattice parameter a and 
the cell angle a seemingly progress together. Thus the 
cell reaches its high-symmetry shape "before" the tran- 
sition to sc takes place; since both a and z must reach 
their high-symmetry values for the transition to occur, 
the A7 — > sc transition pressure is consequently deter- 
mined by the atomic positional parameter z. Thus we 
see that z reaches 0.25 at approximately 21 GPa in the 
LDA and at approximately 28 GPa in the GGA-PBE. 



We consider next, for the LDA and the GGA-PBE, 
the energy- volume and pressure- volume curves of the A7 
and sc structures over a pressure range of approximately 
0-45 GPa, which includes only the A7 — > sc phase tran- 
sition (Fig. |4]). In all cases, the A7 and sc curves merge 
once the A7 — > sc phase transition has occurred, in other 
words, once the volume of the cell has been compressed 
beyond a particular value. This contradicts the recent 
results of Durandurdui^ However, the sc — > A7 trans- 
formation is the result of two continuous distortions, so 
we would expect this transition also to be smooth and 
continuous. We do not expect to see metastable states, 
and we observe none — this agrees with the experi- 
mental findings both of Beister, et air- and of Kikegawa 
and Iwasaki^ yet it disagrees with the theoretical find- 
ings of Da Silva, et ali^ We observe no discontinuities in 
these curves and so our results suggest that the A7 sc 
phase transition is of second order. Experimentally, both 
Beister, et al^ and Kikegawa and Iwasaki^ report that 
the A7 — > sc transition is of first order, where to within an 
experimental error of 0.5% no volume discontinuity is de- 
tected in the former study, whereas a volume discontinu- 
ity of 5% across the transition is found in the latter study. 
Da Silva, et alJ^ report a transition that is weakly first 
order, and a volume discontinuity of less than 1%. We 
believe that if this transition appeared in the past to be 
of first order it was because of the coarse grained nature 
of the investigations performed at the time. 

On the matter of the discrepancies that exist be- 
tween our study and, for example, those of Duran- 
durdu^S and Da Silva, et al.^ it is clear that in those 
investigations the fc-point sampling of the Brillouin zone 
was insufficient: the most dense grid employed by 
Da Silva, et. a/ji^intheirstudy of arsenic was 13x13x13 
— Durandurdu^Si uses only the gamma point for a unit 
cell containing 250 atoms, roughly corresponding to using 
a 5 X 5 X 5 A;-point grid for for a two-atom unit cell. Our in- 
vestigations reveal that these calculations could not have 
been converged. Moreover, in addition to ensuring that 
the Brillouin zone is properly sampled, we believe that 
choosing an appropriate amount of smearing is also very 
important (once again consult the supplementary mate- 
rial for an indepth study of the convergence properties 
with respect to /c-point sampling and smearing of the 
A7 sc transition of arsenic) 

Another point to be made about the insets of Fig. [51 
is that at pressures just below the A7 — > sc transition, 
the nearest neighbor distance actually increases with 
pressure. This occurs over the range of approximately 
15-20 GPa in the LDA and approximately 22-27 GPa in 
the GGA-PBE. The onset of this feature coincides with 
the point at which the a, a and z versus pressure curves 
all begin to experience a kink or slight change of direction 
(see left-hand panels of Fig. [3]) . These pressure intervals 
also correspond to the regions where the A7 pressure- 
volume curves of Fig. 3] start dipping down (from right 
to left) toward those of the simple cubic phase. As the 
volume does not increase discontinuously across the tran- 
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sition, there must be significant atomic restructuring go- 
ing on for the nearest neighbor distance to increase as 
the transition is approached. Although their nearest 
and next-nearest neighbor distances were estimates cal- 
culated according to an empirical relation for the atomic 
positional parameter, Beister, et al^ noticed this point 
as well and noted that such behavior is not unexpected, 
as there is a tendency in covalently bonded solids to dis- 
play increased nearest neighbor distances when undergo- 
ing pressure-induced breaking of directional bonds, ac- 
companied by increased coordination numbers. 

We have performed fits of the energy- volume curves us- 
ing both the third-order Birch-Murnaghan^ and Vinel4£ 
equations to yield values for the bulk modulus. Bo, pres- 
sure derivative of the bulk modulus, B'^, and equilibrium 
volume, Vo, of arsenic for both the A7 and sc structures. 
We find that for the A7 phase, the values we obtain for Bo 
and B'^ are extremely sensitive to the number of points 
included in the fit; this was also found to be the case for 
Needs, et al^ It is also somewhat the case for the sc por- 
tion of the A7 curve, so we have quoted the values that 
we obtain when we fit curves resulting from compressions 
of the exact sc structure, rather than from geometry op- 
timizations. Table HIl summarizes our findings, including 
the results we obtained for the equilibrium cell parame- 
ters from our geometry optimizations of arsenic at zero 
pressure, and compares them to those published in pre- 
vious studies. 



V. CONCLUSIONS 

We have performed simulations of the two-atom unit 
cell of arsenic under compression to investigate its high- 



pressure behavior and compare with experiment. In the 
matter of the A7 — > sc transition, our results strongly 
support the experimental findings of Beister— over those 
of Kikegawa and Iwasaki;^ we furthermore find this tran- 
sition to be of second order. 

Our main critique of the current literature is that we 
believe that any results published until now concern- 
ing this scmi-mctal to metal phase transition arc un- 
convergcd (as no thorough convergence studies have ever 
been published), and that any agreement with experi- 
ment has merely been fortuitous. This study has enabled 
us to present reliably converged results for the A7 sc 
transition of arsenic. 

Using DFT to study any pressure-induced structural 
phase transition involving a metal demands convergence 
testing similar to that which we have carried out^ for 
such purposes, fc-point sampling and smearing must be 
considered thoroughly. 
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TABLE II: Lattice parameters for the uncompressed arsenic structure (A7), in addition to the equihbrium volume (Vo), bulk 
modulus (Bo) and pressure derivative of the bulk modulus (B'o). Those for the uncompressed simple cubic structure of arsenic 
are listed further below. Parametrizations of the generalized gradient approximation to the exchange-correlation functional that 
appear below are PBE,— PW86^ and PB.— (Abbreviations: PP - pseudopotential approach, LAPW - linear augmented 
plane wave approach, BM - fitted with the third order Birch-Murnaghan^ equation, V - fitted with the Vinel4^ equation.) 



A7 structure 


a (A) 


a n 


z 


Vo (A3) 


B, (GPa) 


B'o 


Theory: 














This work, LDA 


3.902 


56.89 


0.231 


19.48(1)™'^ 


60.9(1.6)™'^ 


4.6(4)'^'^, 4.8(4)^ 


This work, PBE 


4.189 


53.96 


0.227 


22.22(5)'^'^'^ 


36.9(1.6)™'^ 


6.1(7)™, 6.25(55)^ 
6.2(5)'^^, 6.3(4)^ 


This work, PW91 


4.195 


53.90 


0.227 


22.30(2)'^'^'^ 


36.9(1.6)'''^'^ 


Ref. W (PP, LDA) 


4.017 


56.28 


0.230 


20.95 


43 




Ref. 44 (PP, LDA) 






0.24 








Ref. 11 (LAPW, LDA) 


4.084 


55.9 


0.2294 


21.8 


77 




Ref. 37 (PP, LDA) 


4.0027 


56.24 


0.230 


20.70 


52 




(PP, PB) 


4.32035 


52.44 


0.225 


23.45 


36 








oz.Ul 


n 00/1 

U.ZZ4 








Ref. (PP, LDA) 


4.031 


56.3 


0.230 


21.18 


59(1) 


4.2(3) 


Ref. 20 (PP, LDA) 








23.57 


64.49 


3.99 


Ref. 21 (PP, PBE) 


4.182 


53.13 


0.225 


21.70 


38.4 


4.34 


Ref. 22 (LAPW, PBE) 


4.246 


53.31 


0.226 


22.82 






Experiment: 














Ref. 34 


4.1018(1) 


55.554(i; 


1 0.22764(4) 


21.303 






Ref. 45 






0.22728(7) 








Ref._3 


4.133 


54.12 


0.227 


21.53 


55.6 


4.4 


Ref. 2 




54.13 






58(4) 


3.3(4) 


Ref.^ 










56(3) 


3.7(2) 


Simple cubic structure 
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Bo (GPa 


) 


B'o 


Theory: 














This work, LDA 


2.635 




^g3QBM,V 


94.53(16)^'^, 93.57(32)^ 


4.29(2)=^^, 4.48(3)^ 


This work, PBE 


2.719 




2o^^qBM,v 


78.99(3)^^^, 78.51(16)^ 


4.318(4)^^, 4.55(2)^ 


This work, PW91 


2.722 




20.17^^'^ 


78.84(4)^^, 78.32(19)^ 


4.297(3)'^'^, 4.53(2)^ 


Ref. (PP, LDA) 


2.687 




19.4 


104 




4.4 


Ref. 44 (PP, LDA) 


2.68 




19.25 


122 




2.32 


Ref. 11 (LAPW, LDA) 


2.717 




20.06 


87 






Ref. (PP, LDA) 


2.704 




19.77 


92(3) 




4.2(3) 


Ref. 20 (PP, LDA) 


2.822 




22.47 


74.15 




3.2 


Ref. 21 (PP, PBE) 


2.674 




19.12 


88.2 




3.93 


Ref. 22 (LAPW, PBE) 


2.731 




20.37 


78.107 




4.317 



